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Abstract 

We investigate the symmetry breaking role of noise in adiabatic quantum computing using the 
example of the CNOT gate. In particular, we analyse situations where the choice of initial Hamil- 
tonian produces symmetries in the Hamiltonian and degeneracies in the spectrum. We show that, 
in these situations, the conventional stipulation that the initial and problem Hamiltonians do not 



> 

commute is unnecessary as noise will inherently play the role of a universal symmetry breaking 



perturbation and split any level crossings that may impede or obstruct the computation. The 
effects of an artificial noise source with a tailored time-dependent amplitude are also explored and 



i-C it is found that such a scheme could offer a considerable performance enhancement. These results 

> are found using a novel, generalised version of the Pechukas- Yukawa model of eigenvalue dynamics. 
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I. INTRODUCTION 



Recently there has been a lot of interest in alternative paradigms to the standard approach 
to quantum computing (i.e. the quantum circuit model). Adiabatic quantum computing 
(AQC) is a promising example which is particularly suited to solving optimisation problems 
PP. AQC involves slow adiabatic evolution from a configuration with an easily reachable 
ground state to one where the ground state encodes the solution to the hard computational 
problem in hand. This scheme is believed to have a number of advantages over the "stan- 
dard" approach, namely, the precise time- dependent control of individual qubits is no longer 
required, and it benefits from an inherent robustness against some of the effects of decoher- 
ence by remaining in the instantaneous ground state at all times [2H1]. Crucially, AQC has 
been shown to be polynomially equivalent, under certain conditions, to the standard gate 
model of quantum computing [6J [T] . The effects of noise on AQC are generally considered to 
be detrimental but manageable [2HI]. ^ was nevertheless stated that its effect can increase 
the success probability of AQC in some situations by; providing an alternative evolution 
trajectory [2], or by thermal relaxation back to the ground state [5J. Here, we investigate 
a more general effect of noise on AQC. Namely, how noise inherently breaks any hidden 
symmetries in the Hamiltonian and, thus, splits any level crossings in the energy spectrum 
which could impede or even prevent the computation. 

AQC can be described by the following Hamiltonian; 

U(\(t)) = H + X(t)ZH bl (1) 

where the ground state of the final Hamiltonian, H , encodes the desired solution, ZHf, is a 
large bias term (with Z ^> 1) and the initial configuration, "H(A = 1) = Hq + ZHj,, has an 
easily reachable, non-degenerate ground state. In order for there to be a high probability 
of the system remaining in the ground state as the bias is switched off, the rate of change 
of the control parameter X(t) must be sufficiently slow to suppress excitation via Landau- 
Zener-Stiickelberg tunneling [5HTU]. 

In AQC, the initial and final Hamiltonians are usually chosen such that no symmetries 
exist in "H(A(t)) by ensuring that they do not commute. This is done to ensure that there 
are no degeneracies in the energy spectrum during the evolution. However, this restriction 
on the choice of initial Hamiltonian may prove practically impossible to realise for generic 
problems. In [11] the example of the adiabatic equivalent of the CNOT gate is discussed as 
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an example of a prototypical quantum algorithm. It was shown that, the choice of a generic 
H b could lead to an abundance of level crossings in the energy spectrum for the 1 00) — > 1 00) 
operation of the CNOT gate. The influence of these level-crossings on the performance of 
the adiabatic quantum algorithm was not considered there. However, the effects of level 
crossings in the spectrum were discussed in p] and it was noted that the addition of an 
appropriate sort of perturbation to the system will break the symmetries and therefore split 
any level crossings. The question of how to provide an appropriate perturbation still remains 
open. 

In this paper, we propose that noise will inherently fulfil the role of the crucial symmetry 
breaking perturbation in physical implementations of an AQC system and as a result of 
this, the condition that H and H b do not commute is no longer required. The CNOT gate 
algorithm is again used as an example of a prototypical quantum algorithm and we show 
that the performance of this generic algorithm is relatively resistant to the effects of noise, 
in agreement with [2HI]. As reported in [2], we find that in some situations the presence of 
noise increases the success probability, we then go on to explore the relationship between 
this increase in success probability and the fidelity of the final state. We also discuss the 
idea of using a tailored artificial noise signal to try and enhance the performance of the 
adiabatic quantum computation process. To do this, we derive and utilise a generalised 
stochastic version of the Pechukas- Yukawa equations p21 [32]; where the dynamics of the 
energy eigenvalues are mapped exactly on to the classical dynamics of a ID gas of Brownian 
particles with a mutual repulsive force. 



II. GENERALISED PECHUKAS- YUKAWA MODEL 

The standard Pechukas- Yukawa model is derived from a Hamiltonian of the form ([I]). 
However, to incorporate a source of noise in to the model we start with the following Hamil- 
tonian; 

H(\(t)) = H + \(t)ZH b + 6h(\{t)), (2) 

where the perturbation strength X(t) plays the role of 'time' and the new stochastic term 
5h(X) describes random fluctuations in the Hamiltonian due to an external noise source. The 
instantaneous eigenvalues and eigenfunctions of (|2| are denoted x n (X) and |n(A)) respec- 
tively; H(X) \n(X)) = x n (X) \n(X)). By following the same procedure as the derivation of the 
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standard Pechukas- Yukawa model, as detailed in [14J, we arrive at the following generalised 
system of equations; 
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where v n (X) = (n\ ZHf, \n) and l nm = (x n — x m ) (n\ ZH^ \m) for n ^ m. These generalised 
equations can be used to describe a wider array of physical systems because of the inclusion 
of the 5h(X) term without any additional assumptions or approximations. Note that if 
the noise term Sh(X(t)) is identically zero the system of equations ^ simply reduces to 
the standard Pechukas- Yukawa equations. The equations ^ describe the dynamics of the 
energy eigenvalues of the Hamiltonian ([2]), but also correspond to the classical dynamics of 
a ID interacting gas, where the nth particle has position x n (X) and velocity v n (X) and the 
strength of the inter-particle force between the nth and the mth particles is described by 

'mm 

(A). 

In order to close the system of equations ([3]), we need to consider the nature of the noise 
term Sh(X). In general, noise in any physical system arises from a number of independent 
sources and therefore as a consequence of the central-limit theorem it seems reasonable to 
assume that the sum of their effects will result in a random Hamiltonian with independent 
Gaussian distributed elements. Such a Hamiltonian will be drawn from one of the Gaussian 
ensembles of random matrix theory. Recent measurements of the low frequency flux noise 
in superconducting flux qubits exhibit a coloured noise spectrum [15]. We therefore assume 
that the noise term, 5h(X), evolves in time as a Ornstein-Uhlenbeck type process, which is 
a simple example of a random process with a coloured spectrum: 

Sh{X) = -r5h(X) + er)(X) (4) 

where r is a correlation time, e is the noise amplitude, and 77(A) is a random matrix valued 
stochastic process where (i](X)) = and (^(A^A')) = <5(A — A'). 



III. THE CNOT GATE 



As mentioned previously AQC has been shown to be polynomially equivalent to the circuit 
model of quantum computing.The 'ground state quantum computing' (GSQC) formalism, 
described in [7], offers the most practical method of constructing a Ho that encodes an 
arbitrary M qubit, N step quantum circuit. In the GSQC formalism, each of the M qubits 
in the circuit is viewed as a single electron that can occupy the states in an array of 2 x (A+l) 
quantum dots; where the rows in the array represents either the |0) or |1) states of the qubit. 
The state of the mth qubit during the nth step of the algorithm is given by the probability 
amplitude of the electron being found on the quantum dots denoted by the indices (m, n, 0) 
and (m,n, 1). This theoretical construction only incurs a polynomial overheard (O(N)) in 
hardware. 

The CNOT gate is one of the simplest entangling quantum gates and when used in con- 
junction with single qubit rotations it forms a set of universal gates. Therefore, any quantum 
algorithm can be constructed using a combination of these gates, because of this property 
we assume that the adiabatic CNOT gate is a representative example of a prototypical adi- 
abatic quantum algorithm. To construct Hq for a CNOT gate using the GSQC formalism, 
we envisage an array of 8 quantum dots (as shown in the inset of Fig. [3]), which corresponds 
to a system of 4 physical qubits. Following the procedure described in |7J we can then write; 

Hcnot = (^ c oioCii ~~ c oooCio) (Cn c oio — CioCooo) 

+ ( c 011^11 — Cooi^lO ^) (Cll c 011 — CT x CioCooi) 

+ CIqCooC^Cu + C01C01C10C10, (5) 

where (^ mn j is a fermionic creation operator that creates an electron on the corresponding 
quantum dot, C^ n = {^ mn0 , ^ mn i\ and a x is a Pauli matrix. The ground state energy of the 
CNOT Hamiltonian Q is zero. To specify the initial state of the qubits before the gate op- 
eration we add a small energy penalty to ([5]); e.g. of the form H Init = u ^c 000 c oo + c Ioo c ioo) 
for the operation 1 00) — > 1 00) . 

We numerically solve the equations ^ and Q. For stability reasons we use the multi- 
step Adams-Moulton method to solve ([3]). For the sake of generality, the initial conditions 
for the Pechukas gas, x n (X = 1), v n (\ = 1) and Z nm (A = 1), are calculated using a per- 
turbation theory expansion in terms of Z~ x as at A = 1, = H + ZH b + 5h(l) = 



Z (Hb + Z~ X H.§ + Z~ l 8h{\)). We assume that has a non-degenerate, well spaced en- 

ergy spectrum. Note that the initial conditions contain all the information about the final 
Hamiltonian H . The initial noise term Sh(X = 1) will be a random matrix drawn from the 
GOE with amplitude e. Throughout the paper, we take /x = —0.1, Z~ x = 0.1 and r = 0.1. 

IV. ENERGY SPECTRA 

Figure []Ja) shows the energy spectra of the CNOT gate acting on the 1 00) and |11) com- 
putational basis states. In both cases, the results from the eigenvalue dynamics simulations 
agree with the results of direct diagonalisation of H Q to four significant figures. In both of 
the spectra, there is an abundance of level crossings which arise because of the symmetries 
of H(\{t)). 

Degeneracies between the ground and first excited states occur in the |11) — > 1 10) op- 
eration, which will result in a success probability of zero for this ideal case. This case can 
be viewed as an example of AQC with an "improper" choice of initial configuration, i.e. 
where [H ,Hb] = 0. It is clear that the addition of any type of perturbation (e.g. noise) will 
break these symmetries, therefore splitting the degeneracies and resulting in a finite success 
probability. The results of the addition of noise are shown in plot Fig. [l](b) and it is evident 
that there are now no degeneracies in either of the spectra. The spectra for the 1 01) — > 1 01 ) 
and 1 10) — > 1 11) operations were found to show similar trends to the |11) — > 1 10) , in that 
they also exhibited a degeneracy between the ground and first excited states in the absence 
of noise. 

V. THE EFFECTS OF NOISE ON LEVEL OCCUPATION STATISTICS 

To properly characterise the effects of noise on the adiabatic quantum algorithm for the 
CNOT gate it is necessary to look at its effects on the probability of level occupation and in 
particular the success probability of the algorithm (i.e. P(n = 0; X(t) = 0\n = 0; A(t) = 1)). 
During the evolution of an adiabatic quantum computer the main mechanism by which the 
system can tunnel from one state to another is Landau-Zener-Stiickelberg tunneling [8HT0]. 
This occurs when the separation of two adjacent energy levels is at a local minimum (i.e. at 
an avoided or level crossing) and the probability of excitation from \n) to \n + 1) via this 
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FIG. 1. (a) The energy spectra of |00) ->■ |00) and |11) ->■ |10) operations of the adiabatic CNOT 
gate in the absence of noise. The solid grey diamonds show the results of direct diagonalisation of 
the Hq. (b) Examples of the spectra for |00) — > |00) and |11) — > 1 10) operations of the CNOT gate 
with a GOE random matrix noise term in the Hamiltonian with amplitude e = 0.1. For illustrative 
purposes, different line styles are used to denote the different energy levels and the third plot in 
each of the subfigures shows an enlarged view of a region of the |11) —> 1 10) spectrum containing 
a number of crossings. 



mechanism is given by 
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, where A min is the minimum separation between x n and x n+ \. We assume that the system 
undergoes uniform evolution (i.e. |A| = where T is the computation time) and that the 
system can initially be found in the ground state with certainty (i.e. P(n = 0; X(t) = 1) = 1). 
Given this information, it is possible to calculate the level occupations as a function of time 
by identifying all the avoided or level crossings in the spectrum and applying the LZS 
equation, ^ to them. 

Figure [2] shows the dependence of the average success probability on the computation 
speed at a range of noise amplitudes for the 1 00) — >■ 1 00) and |11) — > 1 1 0) operations. In 
general, the success probability scales polynomially with the speed, i.e. P oc T~ 7 with 7 ~ 1, 
before approaching unity asymptotically. The exponent is different to those found in [UJ for 
Hamiltonians drawn from the GUE, however there is no reason to expect a single exponent to 
hold universally for all choices of H . From a computational performance point of view, it is 
important to note that the scaling exponent is independent of the noise amplitude, i.e. only 
a prefactor change. This prefactor change could actually be viewed as being beneficial as we 
can see that the success probability at a given computation speed increases linearly with the 
noise amplitude (also shown in Fig. [3J. This occurs as A m j n increases proportionally with e 
and hence this effect is more pronounced for the |11) — > |10) operation where the existence 
of the ground state gap is purely due to noise-induced level splitting. As the |01) — > 1 01) 
and 1 10) — > 1 11) operations also have degenerate ground state gaps, their success probability 
curves behave in a very similar way to the plot shown in Fig. |2]for the 1 11) — > 1 10) operation. 

The increase in success probability with the noise amplitude (at a specific speed) may not 
be wholly beneficial though, as it will come at the cost of the fidelity of the final state (i.e. 
F = \ (0ideai(X = 0)|0 no j se (A = 0))|, where < F < 1). This is because the noise fluctuations 
will drive the state away from the ideal (i.e. in the absence of noise) path of evolution. The 
dependence of the fidelity on noise amplitude is shown in Fig. [3j The fidelity of the final 
state is an important quantity which needs to be maximised to ensure that readout will 
yield the desired solution. For the three operations where the existence of the ground state 
gap is due to noise-induced level splitting, there will be a noise amplitude that offers the 
optimal compromise between the success probability and the fidelity at a given speed and 
this is shown by the intersections of the curves in Fig. [3] These results suggest that the 
conventional stipulation that [Hq,!!^ ^ is unnecessary in practical realisations, as when 
the effects of a generic noise source are taken in to consideration, any hidden symmetries in 
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FIG. 2. Plot of the success probability, averaged over a number of noise realisations, against 
computation speed for the 1 00) — > 1 00) and |11) — > 1 10) operations of the CNOT gate at various 
noise amplitudes. Fits to the polynomial regions, i.e. where P oc T~ 7 , of the curves yield exponents 
of; 7 = 4/3 for the e = case of the 1 00) — > 1 00) operation, and 7 w 1 for all other cases. 

?{(A(t)) will be broken naturally. 



VI. EFFECTS OF AN ARTIFICIAL NOISE SOURCE 



So far we have viewed 5h(X) as the natural effect of a number of noise sources that are 
coupled to the system. If we now envisage a physical system with a negligible level of intrinsic 
noise. In this situation, it may be beneficial to add an artificial random perturbation to the 
system. This would be in order to break any degeneracies in the spectrum and offer an 
alternate, and possibly more efficient, path between the initial and final Hamiltonians. A 
perturbation term with a time-dependent amplitude, which is large enough to widen the 
energy gaps at avoided crossings throughout the majority of the computation process but 
then tends to as X(t) —> 0, would be preferable. For example we could take 



e(A) = e tanh(aA) 
9 
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FIG. 3. Plot showing the trade-off between the fidelity of the final state (dashed line and crosses) 
and the average success probability (solid lines and shapes) as functions of noise amplitude for the 
1 01) — > 1 01) operation. The inset schematic shows the GSQC representation of the 1 01) — > |01) 
operation of the CNOT gate. 

, where a is a constant determining the rate of decay at \(t) > 0. This would ensure that 
5h(X)ij will not be significantly larger than H Qi j at A(t) > 0, where the bias term is small and 
the levels are densely packed. The results of simulations performed using a time- dependent 
amplitude of the form of ^ are shown in Fig. [4j On average, a given success probability 
can be achieved at a much faster computation speed, with an improvement of over 10 2 in 
some cases. In this idealised situation it is also clear to see that the fidelity of the final 
solution state would be unaffected (i.e. F = 1) by this artifical noise signal as e(A = 0) = 0. 



VII. CONCLUSIONS 

We generalise the Pechukas- Yukawa equations to the stochastic case by including an 
additional noise term in the Hamiltonian. This was used in conjunction with a simple, 
yet generic, noise model, based on random matrices and a coloured stochastic process, to 
investigate the effects of noise on the adiabatic algorithm for the CNOT gate. We found 
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(a) |00> to |00> (b)|11>to|10> 




I" 1 

FIG. 4. Plot of the success probability, averaged over a number of noise realisations, against 
computation speed for the 1 00) — > 1 00) and |11) —> 1 10) operations of the CNOT gate with an 
artificial noise source with a time-dependent amplitude and a = 10. 



that, in general, the success probability of the algorithm scaled polynomially as a function 
of computation speed and this scaling was independent of the amplitude of the noise. We 
demonstrate that when the effects of noise are taken in to account, the criteria used to 
select an initial configuration for the system may be relaxed and it is not necessary to 
avoid symmetries in the Hamiltonian. This is because, the presence of noise will break any 
degeneracies in the energy spectrum, crucially those that exist between the ground and first 
excited states. In these situations, the success probability at a given computation speed was 
found to increase linearly with noise amplitude. It was noted that this increase comes at the 
expense of the fidelity of the final state, but an optimal compromise between the two factors 
exists. The effects of an artificial noise source with a time-dependent noise amplitude were 
also investigated. This scheme was found to offer significantly higher success probabilities 
at relatively fast computation speeds and could be engineered in such a way that the fidelity 
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of the solution state is unaffected. 
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